Nuclear matter symmetry energy and the neutron skin thickness of heavy nuclei 
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Correlations between the thickness of the neutron skin in finite nuclei and the nuclear matter 
symmetry energy are studied in the Skyrme Hartree-Fock model. From the most recent analysis 
of the isospin diffusion data in heavy-ion collisions based on an isospin- and momentum-dependent 
transport model with in-medium nucleon-nucleon cross sections, a value of L — 88 ± 25 MeV for 
the slope of the nuclear symmetry energy at saturation density is extracted, and this imposes 
stringent constraints on both the parameters in the Skyrme effective interactions and the neutron 
skin thickness of heavy nuclei. Predicted thickness of the neutron skin is 0.22 ± 0.04 fm for 208 Pb, 
0.29 ± 0.04 fm for 132 Sn, and 0.22 ± 0.04 fm for 124 Sn. 

PACS numbers: 25.70-z, 21.30.Fe, 21.10.Gv 



I. INTRODUCTION 

The study of the equation of state (EOS) of isospin 
asymmetric nuclear matter, especially the nuclear sym- 
metry energy, is currently an active field of research 
in nuclear physics Q, 0, II H E Although 
the nuclear symmetry energy at normal nuclear matter 
density is known to be around 30 MeV from the em- 
pirical liquid-drop mass formula llOf . its values at 
other densities, especially at supranormal densities, are 
poorly known 0, Q • Advances in radioactive nuclear 
beam facilities provide, however, the possibility to pin 
down the density dependence of the nuclear symmetry 
energy in heav y io n collisions induced by t hese nuclei 

H3, HE Ha H3, Ha, Ha, Ha H3, H3, H^, Indeed, signif- 
icant progress has recently been made in extracting the 
information on the density dependence of nuclear sym- 
metry energy from the isospin diffusion data in heavy- 
ion collisions from the NSCL/MSU [H HI HI- Using 
an isospin- and momentum-dependent IBUU04 transport 
model with in-medium nucleon-nucleon (NN) cross sec- 
tions, the isospin diffusion data were found to be consis- 
tent with a relatively soft nuclear symmetry energy at 
subnormal density |37| . 

Information on the density dependence of the nuclear 
symmetry energy can in principle also be obtained from 
the thickness of the neutron skin in heavy nuclei as 
the latter is strongly correlated with the slope L of the 
nuclear matter symmetry energy at saturation density 
H HI HI US El E3 Because of the large uncertainties 
in measured neutron skin thickness of heavy nuclei, this 
has not been possible. Instead, studies have been carried 
out to use the extracted nuclear symmetry energy from 
the isospin diffusion data to constrain the neutron skin 
thickness of heavy nuclei H3, H3- Using the Hartree- 
Fock approximation with parameters fitted to the phe- 
nomenological EOS that was used in the IBUU04 trans- 



port model to describe the isospin diffusion data from the 
NSCL/MSU, it was found that a neutron skin thickness 
of less than 0.15 fm [HHil for 208 Pb was incompatible 
with the isospin diffusion data. 

In the present work, we study more systematically the 
correlation between the density dependence of the nu- 
clear symmetry energy and the thickness of the neutron 
skin in a number of nuclei within the framework of the 
Skyrme Hartree-Fock model. Extracting the values of L 
from the most recently determined density dependence 
of the nuclear symmetry energy from the isospin diffu- 
sion data in heavy-ion collisions, we obtain stringent con- 
straints on the neutron skin thickness of the nuclei 208 Pb, 
132 Sn, and 124 Sn. The extracted value of L also limits the 
allowed parameter sets for the Skyrme interaction. 



II. NUCLEAR SYMMETRY ENERGY AND 
THE SKYRME INTERACTIONS 

The nuclear symmetry energy E sym (p) at nuclear den- 
sity p can be expanded around the nuclear matter satu- 
ration density po as 



E S ym(p) = -Esym(Po) 
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where L and A sym are the slope and curvature of the 
nuclear symmetry energy at po, i-e., 
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The L and K sym characterize the density dependence 
of the nuclear symmetry energy around normal nuclear 
matter density, and thus provide important information 



2 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 



P (fnf 3 ) 



from NSCL/MSU using the IBUU04 transport model 
with in-medium NN cross sections, these interactions are 
recently shown to give, respectively, the upper and lower 
bounds for the stiffness of the nuclear symmetry energy 
|37| . It is seen from Fig. ^that the density dependence 
of the symmetry energy varies drastically among differ- 
ent interactions. Although the values of E sym (po) are all 
in the range of 26-35 MeV, the values of L and K sym are 
in the range of -50-100 MeV and -700-50 MeV, respec- 
tively. 



III. NEUTRON SKIN THICKNESS OF FINITE 
NUCLEI AND THE SLOPE OF NUCLEAR 
SYMMETRY ENERGY AT SATURATION 
DENSITY 



FIG. 1: (Color online) Density dependence of the nuclear 
symmetry energy E ayln (p) for 21 sets of Skyrme interaction 
parameters. The results from the MDI interaction with x = 
— 1 (open squares) and (solid squares) are also shown. 



on the properties of nuclear symmetry energy at both 
high and low densities. 

In the standard Skyrme Hartree-Fock model [381 HI . 
EE EE EE EH j the interaction is taken to have a zero- 
range, density- and momentum-dependent form and the 
Skyrme interaction parameters are chosen to fit the bind- 
ing energies and charge radii of a large number of nuclei 
in the periodic table. For infinite nuclear matter, the nu- 
clear symmetry energy from the Skyrme interaction can 
be expressed as 0, EE] 



1 h 2 



Es ^ ip) 3 2^ I ~ 
1 



3^ 
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+ 24 ("J [ - 3 ^ 
+ (4 + 5.x 2 )t 2 ]p 5 / 3 , 



(4) 



where the cr, to — 13 , and xo — x 3 are Skyrme interaction 
parameters. 

Fig. n displays the density dependence of E sym (p) 
for 21 sets of Skyrme interaction parameters, i.e., SKM, 
SKM*, RATP, SI, SII, SIII, SIV, SV, SVI, E, E,,, G ff , 
R a , Z, Z a , Z*, T, T3, SkX, SkXce, and SkXm. Values of 
the parameters in these Skyrme interactions can be found 
in Refs. [44l EE EE] ■ F° r comparison, we also show in 
Fig. ^ results from the phenomenological parametriza- 
tion of the momentum-dependent nuclear mean-field po- 
tential based on the Gogny effective interaction EE], i.e., 
the MDI interactions with x = — 1 (open squares) and 
(solid squares), where different x values correspond to 
different density dependence of the nuclear symmetry en- 
ergy but keep other properties of the nuclear EOS the 
same [36]. From comparing the isospin diffusion data 



The neutron skin thickness S of a nucleus is defined as 
the difference between the root-mean-square radii \J (r n ) 
of neutrons and U (r p ) of protons, i.e., 



(5) 



It has been known that S is sensitive to the density de- 
pendence of the nuclear symmetry energy, particularly 
the slo pe p arameter L at the normal nuclear matter den- 
sity UllliEEEEJEl- Usin S above 21 sets of Skyrme 
interaction parameters, we have evaluated the neutron 
skin thickness of several nuclei. In Figs. EJa), (b) and (c), 
we show, respectively, the correlations between the neu- 
tron skin thickness of 208 Pb with L, K sym , and -E S ym(po)- 
It is seen from Fig. Ufa) that there exists an approximate 
linear correlation between S and L. The correlations of S 
with K sym and E sym (po) are less strong and even exhibit 
some irregular behavior. The solid line in Fig. EJa) is a 
linear fit to the correlation between S and L and is given 
by the following expression: 



or 



5( 208 Pb) = (0.1066 ±0.0019) 

+ (0.00133 ±3.76 x 10~ 5 ) x L, 



(-78.5 ±3.2) 

+ (740.4 ±20.9) x 5( 208 Pb), 



(6) 



(7) 



where the units of L and S are MeV and fm, respectively. 
Therefore, if the value for either S'( 208 Pb) or L is known, 
the value for the other can be determined. 

It is of interest to see if there are also correlations be- 
tween the neutron skin thickness of other neutron-rich 
nuclei and the nuclear symmetry energy. Fig. [3] shows 
the same correlations as in Fig. [21 but for the neutron- 
rich nuclei 132 Sn, 124 Sn, and 48 Ca. For the heavy 132 Sn 
and 124 Sn, we obtain a similar conclusion as for 208 Pb, 
namely, S exhibits an approximate linear correlation with 
L but weaker correlations with K sym and -E S ym(/?o)- For 
the lighter 48 Ca, on the other hand, all the correlations 
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TABLE I: Linear correlation coefficients Ci of S with L, K ayul 
and E sym (p ) for 208 Pb, 132 Sn, 124 Sn, and 48 Ca from 21 sets 
of Skyrme interaction parameters. 



Ci (%) 


208 p b 


132 Sn 


124 Sn 


48 Ca 


S-L 


99.25 


98.76 


98.75 


93.66 




92.26 


92.06 


92.22 


86.99 


S-Esym 


87.89 


85.74 


85.77 


81.01 
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FIG. 2: (Color online) Neutron skin thickness S of 208 Pb 
as a function of (a) L, (b) K aym , and (c) E synl (po) for 21 
sets of Skyrme interaction parameters. The line in panel (a) 
represents a linear fit. 



nuclei is determined by the pressure difference between 
neutrons and protons, which is proportional to the pa- 
rameter L. 

To give a quantitative estimate of above discussed cor- 
relations, we define the following linear correlation coef- 
ficient C;: 
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FIG. 3: (Color online) Same as Fig. 2 but for nuclei 132 Sn 
(Solid squares), 124 Sn (Open squares) and 48 Ca (Triangles). 



become weaker than those of heavier nuclei. Therefore, 
the neutron skin thickness of heavy nuclei is better cor- 
related with the density dependence of the nuclear sym- 
metry energy. As in Eq. @J and (7J, a linear fit to the 
correlation between S and L can also be obtained for 
132 Sn and 124 Sn, and the corresponding expressions are 

S( 132 Sn) = (0.1694 ± 0.0025) 

+ (0.0014 ± 5.12 x 10~ 5 ) x L, (8) 

L = (-117.1 ±5.4) 

±(695.1 ± 25.3) x 5( 132 Sn), (9) 



and 



S'( 124 Sn) = (0.1255 ± 0.0020) 

±(0.0011 ± 4.05 x 10~ 5 ) x L. (10) 

L = (-110.1 ±5.2) 

±(882.6 ± 32.3) x 5( 124 Sn), (11) 



Similar linear relations between S and L are also ex- 
pected for other heavy nuclei. This is not surp rising as 
detailed discussions in Refs. 0, H HI ES El El have 
shown that the thickness of the neutron skin in heavy 



In the above, A and B are the linear regression coeffi- 
cients, (xi, yi) are the sample points, and n is the num- 
ber of sample points. The linear correlation coefficient 
Ci measures the degree of linear correlation, and C/ = 1 
corresponds to an ideal linear correlation. Table [I] gives 
the linear correlation coefficient Ci for the correlation of 
S with L, K sym and E sym ( Po ) for 208 Pb, 132 Sn, 124 Sn, 
and 48 Ca shown in Figs. |5] and [3] for different Skyrme 
interactions. It is seen that these correlations become 
weaker with decreasing nucleus mass, and a strong lin- 
ear correlation only exists between the S and L for the 
heavier nuclei 208 Pb, 132 Sn, and 124 Sn. Therefore, the 
neutron skin thickness of these nuclei can be extracted 
once the slope parameter L of the nuclear symmetry en- 
ergy at saturation density is known. 



IV. CONSTRAINTS FROM ISOSPIN 
DIFFUSION DATA IN HEAVY ION COLLISIONS 

Experimentally, the degree of isospin diffusion between 
the projectile nucleus A and the target nucleus B can be 
studied via the quantity Ri [35l IHfij , 



R,. = 



2X 



A+B 



-X A + J 



X 



B+B 



X B+B 



(15) 



where X is any isospin-sensitive observable. By con- 
struction, the value of Ri is 1 (—1) for symmetric 
A + A (B + B) reaction. If isospin equilibrium is reached 
during the collision as a result of isospin diffusion, the 
value of Ri is about zero. In the NSCL/MSU exper- 
iments with A = 124 Sn and B = 112 Sn at a beam 
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energy of 50 MeV/nucleon and an impact parameter 
about 6 fm, the isospin asymmetry of the projectile- like 
residue was used as the isospin tracer X j^J. Using 
an isospin- and momentum-dependent IBUU04 transport 
model with free-space experimental NN cross sections, 
the dependence of Ri on the nuclear symmetry energy 
was studied from the average isospin asymmetry of the 
projectile-like residue that was calculated from nucle- 
ons with local densities higher than po/20 and velocities 
larger than 1/2 the beam velocity in the center-of-mass 
frame 36]. Comparing the theoretical results with the 
experimental data has allowed us to extract a nuclear 
symmetry energy of E sym (p) « 31.6(p/po) 105 - Includ- 
ing also medium-dependent NN cross sections, which are 
important for isospin-dependent observables |37tl5lj ]. the 
isospin diffusion data leads to an even softer nuclear sym- 
metry energy of E sym (p) ss 31.6(p//0n) 7 with 7 « 0.7 |37| . 

In Fig. 01 we show the results from the IBUU04 trans- 
port model with in-medium NN cross sections, that are 
consistent with the mean-field potential obtained with 
the MDI interactions used in the model, for the degree of 
the isospin diffusion 1 — Ri as a function of L. The shaded 
band in Fig. H indicates the data from NSCL/MSU 
It is seen that the strength of isospin diffusion 1 — Ri 
decreases monotonically with decreasing value of x or in- 
creasing value of L. This is expected as the parameter 
L reflects the difference in the pressures on neutrons and 
protons. From comparison of the theoretical results with 
the data, we can clearly exclude the MDI interaction with 
x = 1 and x — — 2 as they give either too large or too 
small a value for 1 — Ri compared to that of data. The 
range of x or L values that give values of 1 — Ri falling 
within the band of experimental values could in princi- 
ple be determined in our model by detailed calculations. 
Instead, we determine this schematically by using the re- 
sults from the four x values. For the centroid value of L, 
it is obtained from the interception of the line connecting 
the theoretical results at x = — 1 and with the central 
value of 1 - Ri data in Fig. H i.e., L = 88 MeV. The 
upper limit of L = 113 MeV is estimated from the inter- 
ception of the line connecting the upper error bars of the 
theoretical results at x = —1 and —2 with the lower limit 
of the data band of 1 — R4. Similarly, the lower limit of 
L — 65 MeV is estimated from the interception of the line 
connecting the lower error bars of the theoretical results 
at x — and —1 with the upper limit of the data band of 
1 — R4. This leads to an extracted value of L = 88 ± 25 
MeV as shown by the solid square with error bar in Fig. 

El 

The extracted value of L — 88 ± 25 MeV gives a rather 
stringent constraint on the density dependence of the nu- 
clear symmetry energy and thus puts strong constraints 
on the nuclear effective interactions as well. For the 
Skyrme effective interactions shown in Fig. ^ for in- 
stance, all of those lie beyond x — and x = — 1 in 
the sub-saturation region are not consistent with the ex- 
tracted value of L. Actually, we note that only 4 sets of 
Skyrme interactions, i.e., SIV, SV, G CT , and R CT , in the 21 



E/A=50 MeV and b=6 fm 
-□- IBUU04 with MDI + a 




60 90 120 
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180 



FIG. 4: (Color online) Degree of the isospin diffusion 1 — Ri as 
a function of L using the MDI interaction with x — —2, —1, 0, 
and 1. The shaded band indicates the data from NSCL/MSU 
[^3 • The solid square with error bar represents L — 88 ± 25 
MeV. 



sets of Skyrme interactions considered here have nuclear 
symmetry energies that are consistent with the extracted 
L value. 

The extracted L value also allows us to determine from 
Eqs. ©, ©, and Q10p. respectively, a neutron skin thick- 
ness of 0.22 ± 0.04 fm for 208 Pb, 0.29 ± 0.04 fm for 132 Sn, 
and 0.22±0.04 fm for 124 Sn. Experimentally, great efforts 
were devoted to measure the thickness of the neutron 
skin in heavy nuclei [521 153| , and a recent review can be 
found in Ref. (54J. The data for the neutron skin thick- 
ness of 208 Pb indicate a large uncertainty, i.e., 0.1-0.28 
fm. Our results for the neutron skin thickness of 208 Pb 
are thus consistent with present data but give a much 
stronger constraint. A large uncertainty is also found ex- 
perimentally in the neutron skin thickness of 124 Sn, i.e., 
its value varies from 0.1 fm to 0.3 fm depending on the ex- 
perimental method. The proposed experiment of parity- 
violating electron scattering from 208 Pb at the Jefferson 
Laboratory is expected to give another independent and 
more accurate measurement of its neutron skin thickness 
(within 0.05 fm), thus providing improved constraints on 
the density dependence of the nuclear symmetry energy 

Most recently, an accurately calibrated relativistic 
parametrization based on the relativistic mean-field the- 
ory has been introduced to study the neutron skin thick- 
ness of finite nuclei [5JJ. This parametrization can de- 
scribe simultaneously the ground state properties of fi- 
nite nuclei and their monopole and dipole resonances. 
Using this parametrization, the authors predicted a neu- 
tron skin thickness of 0.21 fm in 208 Pb, 0.27 fm in 132 Sn, 
and 0.19 fm in 124 Sn [57], E^]. These predictions are in 
surprisingly good agreement with our results constrained 
by the isospin diffusion data in heavy-ion collisions. 
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V. SUMMARY 

In summary, we have studied the correlation between 
the neutron skin thickness of finite nuclei and the nuclear 
symmetry energy within the framework of the Skyrme 
Hartree-Fock model. As in previous studies, we have 
found a strong linear correlation between the neutron 
skin thickness of heavy nuclei and the slope L of the 
nuclear matter symmetry energy at saturation density. 
This correlation provides stringent constraints on both 
the density dependence of the nuclear symmetry en- 
ergy and the thickness of the neutron skin in heavy nu- 
clei. From the most recent analysis of the isospin dif- 
fusion data in heavy-ion collisions using an isospin- and 
momentum-dependent transport model with in-medium 
NN cross sections, the value L = 88 ± 25 MeV has been 
extracted. The relatively constrained value for the slope 
of the nuclear matter symmetry energy imposes strong 



constraints on the parameters in the Skyrme effective in- 
teractions and leads to predicted neutron skin thickness 
of 0.22 ± 0.04 fm for 208 Pb, 0.29 ± 0.04 fm for 132 Sn, and 
0.22 ± 0.04 fm for 124 Sn. 
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